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TRANSPORT  THEORY  FOR  HEAVY  ION  BEAMS  NEAR  A  TARGET  SURFACE 


I.  INTRODUCTION 

When  treated  from  the  point  of  view  of  transport  theory,  the  bombardment 

of  a  target  by  a  beam  of  energetic  particles  involves  highly  singular 

functions,  where  by  singular  functions  we  mean  either  Dirac  delta  functions 

or  functions  that  are  sharply  peaked  in  one  or  more  of  the  independent 

variables.  At  the  penetrating  surface  the  distribution  function  for  the 

penetrating  particles  has  Dirac  delta  functions  in  energy  and  angle 

corresponding  to  the  monoenergetic  and  paraxial  nature  of  the  bombarding 

beam.  In  this  work  we  undertake  a  systematic  treatment  of  these 

singularities.1  We  show  that,  if  one  is  interested  in  a  numerical 

calculation  of  the  distribution  function,  numerical  approximations  to  the 

singular  functions  can  be  avoided  by  systematically  extracting  them  and 

treating  them  exactly.  We  also  show  how  our  approach  naturally  leads  to  an 

iterative  (successive  scattering)  treatment  of  the  distribution  function 

near  the  penetrating  surface  that  could  be  the  basis  of  a  self-contained 

transport  theory  of  near-surface  phenomena  in  ion  bombardment. 

The  transport  theory  for  penetration  of  an  amorphous  target  is 

formulated  in  the  next  section.  Our  approach,  which  is  based  on  the 

Boltzmann  equation,  is  familiar  to  workers  in  neutron  transport 
2-4 

theory.  Specifically,  our  approach  parallels  and  was  suggested  by  that 

5 

of  Manning  and  Padgett  ( MP ) .  In  contrast,  much  of  the  recent  work  on  ion 

6-9 

bombardment  utilizes  a  transport  equation  of  the  type  used  by  Lindhard, 
Scharff,  and  Schiott  (ISS).1^  Our  treatment  of  the  inelastic  collisions 
is  that  of  ISS. 

In  Section  III  we  show  that  the  distribution  function  can  be  split  into 
a  hierarchy  of  terms  of  decreasing  singularity:  a  %  term,  which 
represents  the  beam  as  it  first  enters  the  target,  wherein  it  loses  energy 
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to  inelastic  processes  but  retains  its  unidirectional  character;  terms 

^  ,  ....  -/  ,  which  represent  beam  ions  that  have  undergone  one, 

two,  ...,  n-1  elastic  collisions;  and  a  term,  which  represents  beam 
ions  which  have  undergone  n  or  more  elastic  collisions.  The  explicit 
solutions  for  the  first  two  terms  are  given  in  Secs.  IV  and  V.  The  formal 
solutions  for  the  remaining  fg  terms,  s  =  2, 3,..., n-1,  and  the  equation 
that  must  be  solved  for  are  discussed  in  Sec.  Vl.  The  following  section 
presents  the  results  of  some  sample  numerical  calculations  for 
and  $ ,  using  the  Lindhard  cross  sections.^  Section  V 1 1 1  contains  our 
concluding  remarks  and  a  discussion  of  a  number  of  possible  extensions  of 
this  work. 

The  conceptual  form  of  the  solution  of  the  Boltzmann  equation  that  we 
use  (Sec.  Ill)  is  the  same  as  that  of  Manning  and  Padgett.  The  solution  for 
the  most  singular  term,  the  term,  is  al_.  the  same  as  that  of  MP, 
although  the  derivation  is  different.  Our  version  of  the  W  term  differs 
somewhat  from  that  of  MP;  it  provides  a  complete  description  of  particles 
that  have  undergone  one  elastic  collision  whi le*  suffering  the  effects  of  the 
continuous  inelastic  losses.  The  iteration  method  of  Sec.  Vl  for  solving 
for  successive  %  terms  is  similar  to  that  discussed  by  MP,  but  the  present 
effort  provides  formal  expressions  for  the  result  of  each  iteration. 

The  continuous  energy  loss  model  for  inelastic  collisions  used  by 
ISS10  is  further  discussed  in  Appendix  A.  Our  presentation  is  based  on 
the  work  of  Symon,^  and  follows  that  of  Ref.  5.  Appendix  B  contains  a 
brief  description  of  the  LSS  cross  section,1^1  using  our  notation. 


II.  BOLTZMANN  EQUATION  FOR  BEAM  PENETRATION 

For  convenience,  we  repeat  here  the  formulation  of  the  Boltzmann 

5 

equation  for  beam  penetration  as  found  in  MP.  Consider  a  beam  of  heavy 

ions  of  charge  Z^,  mass  m^,  and  energy  Eg  penetrating  an  amorphous 

target  made  up  of  atoms  of  charge  Z,  and  mass  m„,  with  number 

c  c  2 

density  N.  Following  the  notation  of  Case  and  Zweifel,  we  take  the 
single  particle  distribution  function  to  be  where 
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X(c.t<)  3  the  number  of  particles  in  Sr  about  c  and  Sts  about  if  at 

time  ( 1 ) 

The  function  TCg'+y)  is  defined  as  follows: 

A/  *'  Sr  SvJt*  V'  =  the  probable  number  of 

beam  particles  of  velocity  in  Si/'  about  £  ,  located  in  Sr 

about  t, 

which  undergo  a  collision  in  time  about  "t  such  that  their  final 
velocities  lie  in  A  about  g.  (2) 

Since  X  WC.K,*)  is  the  flux  of  particles  with  velocity^,  we  see 
that  rlr'-ty)  is  a  collision  cross  section.  The  total  cross  section  is 
given  by 

-  Js*  r(jc'fy)  .  (3) 

(Note  that  Case  and  Zweifel  use  the  various  T  quantities  to  represent 

macroscopic  cross  sections,  which  include  the  factor/^.)  By  applying  the 

principle  of  the  conservation  of  particles  to  the  7^  function  over  a  time 

interval  /f  as  c  is  changea  to  c  ttfg  and  ^changes  to  ,  one 

2  5 

finds  ’  the  Boltzmann  transport  equation  for  beam  penetration  to  be 

+  X-ZTC  +  v»rs(«)  %  = 

*  N  jSS  r(y'-*y)  (*) 

where  j </-  represents  the  source  term  for  the  ion  flux. 

Throughout  the  remainder  of  this  paper  we  will  make  the  simplifying 
assumptions  of  rotational  invariance,  translational  invariance,  plane 
symmetry,  and  time  independence.  We  presume  that  the  beam  direction  is 
parallel  to  the  x-axis  and  that  the  beam  enters  the  target  perpendicular  to 
its  face.  We  let  x  measure  the  distance  of  penetration  from  the  surface  and 
let 

S  *  *  (5) 
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represent  the  cosine  of  the  angle  between  the  beam  direction  in  the  target 
and  the  initial  beam  direction.  With  the  above  assumptions,  the  cross 
sections  depend  only  on  E  ,  and  .  We  define  a  new  function 
that  is  related  to  by 

X <ie  /*  4*.  *  ¥*  X(c.Z,  T)  Jo  (6) 

and  that  satisfies  the  transport  equation 

7x  X(f**/*)  *  A/  rt(£)  *  *f(£,jr,yuJ / if 

+  A//SS  r-fcf,  (7) 

Because  of  the  assumed  time  independence,  it  is  simpler  to  work  with  the 
vector  flux,  rather  than  the  number  density  X(e^).  In  the  case  of  plane 
symmetry,  the  vector  flux  is  given  by 

ft £,*,/«)  •  x*  *  •  *  (8) 

and  it  satisfies  the  Boltzmann  equation 

A  &  +  #  r,(£j  ?  * 

+  CwrVy*'  r (£,£',  }•}')  ?(%*,£>').  0) 

We  now  make  the  widely  used  assumption1^  that,  insofar  as  the 
statistics  of  the  transport  equation  are  concerned,  the  atomic  collisions 
can  be  modeled  as  two  separate  types,  inelastic  and  elastic,  and  that  these 
two  types  of  collisions  occur  independently  and  without  interference. 
Inelastic  collisions  are  thought  of  as  leaving  the  electrons  in  the  target 
atom  in  an  excited  state,  but  with  the  center-of-mass  motion  of  the  target 
atom  unperturbed.  Elastic  collisions  are  thought  of  as  transferring 
momentum  to  the  target  atom  as  a  whole,  while  leaving  unchanged  its 
electronic  state. 

Correspond ing  to  this  assumption,  we  can  split  the  cross  section  into 
two  parts.  For  the  inelastic  losses  we  use  the  continuous  energy  loss 
model,11  for  which  the  energy  loss  results  from  infinitely  many 
collisions,  each  of  which  transfers  an  infinitesimal  energy.  The  direction 
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of  motion  of  the  particles  is  unaltered  by  the  inelastic  losses.  The 
inelastic  cross  section,  which  satisfies  the  relation 

sir)  «  cr-e)  j  ,  no) 

is  discussed  in  Appendix  A.  The  function  see')  is  called  the  electronic  or 
inelastic  stopping  power;  it  represents  the  rate  of  energy  loss  to  inelastic 
processes.  The  modifications  that  this  form  of  continuous  loss  mechanism 
generates  in  the  8oltzmann  equation  appear  in  Appendix  A.. 

In  the  case  of  elastic  collisions,  wherein  the  beam  ion  velocity  goes 
from  %  to  conservation  of  energy  and  momentum  yield 

jCM  *  il  (»+<)(£'/£)*  -(*-<)(*/*')*]  on 

as  the  cosine  of  the  scattering  angle  in  the  laboratory  system,  with 

M  =  .  (12) 

The  maximum  energy  transfer  in  an  elastic  collision  of  particles  of  mass 
m-j  and  m2  is  given  by 

Tmar  =  I*?'  =  Hr*,  (/*,  £’  f  (13) 

where  e'  is  the  energy  of  the  impinging  particle.  Consequently,  the  minimum 
energy  for  the  outgoing  particle  is 

emi„  =  e'/jf  =  £'  ( m. -m) * /  (m,  +**) a  #  (14) 

We  now  write  the  elastic  portion  of  the  cross  section  in  the  form 

m  (e',e,  0  *  M"  f(K£) .  05) 

The  quantity  F(e/e)  contains  the  unit  step  function  due  to  the 

kinematical  considerations.  For  convenience,  we  also  define 

06) 
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The  total  elastic  cross  section  is  now  written 

rr(e')  =  Jfyp  Je  & (K a)  .  (i7) 

In  the  light  of  the  discussion  of  the  last  two  paragraphs,  we  obtain  the 
Boltzmann  equation  in  the  form 

STx  -  t/VrC*)  ¥  -  £(  ¥') 

=  ytt  V (£, *,/*)  *  &  J£  xte'  j (is) 

III.  FORM  OF  THE  SOLUTION 

We  assume  that  the  incident  beam  is  perfectly  coll imated*and 
monochromatic.  In  particular,  we  take  the  source  term  in  Eq.  (18)  to  be 

/*  S(/*-0  S(e-fo)  (19) 

where  6  is  the  Dirac  delta  function.  The  source  has  units  of  number  of 
particles  per  unit  area  per  unit  time.  Since  we  are  examing  a  time 
independent  problem,  all  of  our  results  will  be  proportional  to 
henceforth  we  will  assume  unit  incoming  flux,  i.e.,  ^1*/. 

This  section  outlines  the  form  of  the  solution  that  we 

choose.  The  crucial  element  of  this  analysis  is  that  we  give  special 
attention  to  the  behavior  of  near  the  x=o  plane,  so  that  the  singular 

portions  of  {?  can  be  dealt  with  separately  and  specifically.  Then,  when 
the  remainder  of  $  is  evaluated,  the  strongly  peaked  portions  will  not 
disrupt  the  numerical  processes.  We  consider  here  only  the  distribution 
function  for  the  penetrating  beam.  The  extension  of  this  treatment  to 
include  recoil  contributions  ought  to  be  straightforward.”  Our  approach 
to  this  problem  appears  somewhat  different  from  that  of  Manning  and  Padgett, 
but  conceptually  they  are  the  same  and  our  treatment  evolved  from  theirs. 
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For  the  sake  of  convenience,  we  follow  Winterbon,  Sigmund,  and  Sanders 
(WSS)®  in  presuming  that  the  target  material  is  present  on  both  sides  of 
the  source  plane.  For  cases  where  this  is  a  bad  approximation  to  a  beam  of 
particles  incident  on  a  target  face,  the  corresponding  will  indicate 

considerable  scattering  back  and  forth  across  the  x=o  plane. 

One  way  to  solve  for  We,*,#)  is  to  resolve  it  into  components  according 
to  the  number  of  elastic  collisions  that  the  beam  particles  have  undergone; 
that  is,  we  write 


*  £ 


% 


* 


i. 


(20) 


The  term  represents  incident  particles  that  have  undergone  no  elastic 
collisions,  the  %  terms,  s  =  l,2,...,n  -  1,  represents  those  that  have 
undergone  exactly  s  elastic  collisions,  and  the  term  represents  particles 
that  have  undergone  n  and  more  collisions.  The  decision  as  to  how  many 
specific  %  terms  to  evaluate,  before  attempting  to  solve  for  ,  depends 
on  the  convenience  of  calculating  the  %  quantitites  and  the  degree  to 
which  the  singular  portions  of  have  been  removed. 

If  we  substitute  the  form  (20)  into  Eq.  (18),  with  the  source  (19),  then 
the  Boltzmann  equation  can  be  resolved  into  a  series  of  successive 
scattering  equations.  For  convenience,  we  define  the  operators 

£  s  /*  3x  -  S(S)  (2l) 

Xz  5  Hrr(e)  -  (22) 

so  that  Eq.  (18)  is  equivalent  to 

(l,* Li)  %  (£,*,#)  ~  (£,*y#)  /  3*0,1,...,*-/ ,  (23) 

(L*li)  s  (fS/#)  (24) 


where  , 

*  Jv  J€  Je'  &(?,£) J Jr ‘fUfj***)  (25) 

I/,  (26) 
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ana 


W-t  ( £,*,/<)  -  (f,  *//«)  .  (27) 

For  the  case  of  no  elastic  collisions,  that  is,  F(t:*)-o  ,  all  of  the 
terms  in  i?  will  remain  zero,  except  for  the  %  term.  The  continuous 
energy  loss  model  for  inelastic  collisions  results  in  the  beam  remaining 
collimated  and  monochromatic  as  it  penetrates,  but  with  the  energy  per 
particle  continuously  decreasing.  If  we  now  imagine  that  the  elastic 
scattering  is  turned  on,  when  a  particle  undergoes  its  first  elastic 
collision  it  leaves  the  'K  term  and  enters  the  ft  term.  Eventually,  after 
its  nth  elastic  collision,  it  enters  the  term,  which  represents  it 
thereafter.  Thus,  as  the  beam  penetrates  deeper  into  the  target,  we  expect 
the  terms  to  first  grow  in  magnitude  and  then  attenuate, 

leaving  only  the  J*  term  deep  in  the  target  material. 

The  approach  just  described  is  equivalent  to  the  use  of  a  hierarchy  of 
limiting  procedures  by  Manning  and  Padgett. 

A  comment  on  the  handling  of  the  inelastic  losses  is  in  order  before  we 
proceed  to  the  derivations.  The  results  of  the  next  several  sections  are 
obtained  for  a  general  inelastic  loss  law  S(0) ,  but,  in  addition  to  the 
general  results,  we  present  the  specific  results  for  the  velocity  stopping 
1  aw 


S(£J-  K  F*  ,  (28) 

A  velocity  stopping  law  form  is  given  as  a  primed  equation  following  the 
general  result,  wherever  appropriate.  We  do  this,  firstly,  because  the 
velocity  stopping  law  is  the  form  most  commonly  used  to  represent  inelastic 
losses  at  low  energies  and  it  is  convenient  to  have  the  results  at  hand. 
Further,  the  example  of  this  simple  form  of  S(£)  serves  to  elucidate  the 
method  we  use  of  solve  the  resolved  Boltzmann  equations  (Eqs.  (23-26)). 
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IV.  THE  SOLUTION  FOR  %  (£,*,/*) 


We  will  solve  for  ft  (*,'*")  indeed,  for  each  %(?/*//#)  by  using  the 

12-14 

method  of  characteristics  of  partial  differential  equation  theory. 

We  have  adopted  this  method  with  little  explanation  of  its  use,  so  the 
reader  may  want  to  examine  Refs.  (12-14),  or  others,  before  proceeding. 

The  equation  satisfied  by  is 

*  A  f.  (29) 


The  method  of  characteristics  demonstrates  that  a  first  order  partial 
differential  equation  is  equivalent  to  a  coupled  set  of  ordinary 
differential  equations;  specifically,  Eq.  (29)  is  equivalent  to 


=  JlL  -  -  ^  ,  *4 ft _ 

o  ^  S(SJ  -XiK+J- 


(30) 


The  first  expression  merely  indicates  that  we  may  consider  the  direction 
cosine /e  as  a  parameter  instead  of  as  a  variable.  The  next  pair  of  equated 
terms  in  the  string  (30)  are  easily  solved,  yielding 


JL 


-  Jt  s"(f)  =  A  *  £  -  fie) 


(31) 
(31 ' ) 


which  defines  fie) ,  where  A  is  an  arbitrary  constant  and  £*  is  arbitrary 
except  that  we  take  it  to  be  greater  than  the  beam  energy  £9  .  Clearly, 
f  (£)  will  be  a  known  function,  numerically,  at  least,  since  see)  is  known. 

The  method  of  characteristics  has  thus  indicated  that  if  is  a 
solution  of  Eq.  (29),  then  so  is  j ,  where  S(A)  is  any  arbitrary 
function  of  A.  That  this  is  true  can  be  seen  by  substituting  an  arbitrary 
function  of  »» a-fif)  into  Eq.  (29),  and  thus  verifying  that  it  is  a 
solution  of  the  homogeneous  portion  of  that  equation.  This  fact  allows  us 
to  use  Eq.  (31)  to  eliminate  one  of  the  two  variables  ( x  and  £T)  in  favor  of 
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the  other  when  we  solve  for  %  . 

The  last  equated  pair  of  expressions  in  Eq.  (30)  can  be  written 

=  s"(e)  [  /,  £  -vj  . 

We  introduce  the  integrating  factor  ,  with 

I(t:e)  =  <*p  l  -Je  Jc  S'?*)  S*(t)  } 


so  that 


*  xe’,ej 

»  -  ne’.f)  s"fe)  je  . 


A  I 

If  we  insert  the  explicit  form  (22)  for  J f.a  ,  we  see  that  Z(E/£)  can  be 


written 


X($*,e)  -  see0)  £(?*,£)  s"(e)  (36) 

R(Z%5)  -  tjep  {-J*S P  $"{*)  Mt}(£)  ]  .  (37) 

Upon  using  expression  (19)  for  and  using  Eq.  (31)  to  eliminate  x  in 
favor  of  £  in  Eq.  (35),  we  find 

if  -  B  +jf  Sc  JUJ’)S-(t)  $!»-!) Ucfa)  (38) 

where  3  is  an  arbitrary  constant;  upon  performing  the  integration  we  obtain 

%  *  8  +  U(c,  -  E)  X(ef/  6  *)  S"(fa)  fr-o  i(A  -ftoJ)  .  ( 39 ) 

Our  boundary  conditions  are  simply  that  the  only  particles  of  type  "2"  in 
the  system  are  those  that  arise  from  the  source  «/.  Clearly,  then,  we  have 
that  the  inhomogeneous  contribution  vanishes;  that  is,  we  have  E~0.  By 
using  expression  (31)  we  have,  with  /*•!  ,  that 
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-  x  -  }f*  A  s-CtJ 

(40) 

A* fife)  •  x  -  a(e,*  -S*)/K  . 

(40  * ) 

Let  us  define  &(*)  to  be  a  function  that,  when  substituted  for  £  in 
Eq.  (40),  makes  the  expression  vanish.  Thus  we  can  write,  for  £*£(*)  ^ 

x-J?  Jt  s-(t)  -  (£'*,(*))  s"(euj) . 


(For  the  velocity  stopping  law,  we  have 

e,(x)  -  (£**  -  -k/fx)*  (41 1 ) 

as  the  explicit  form  of  £.(x)  .) 

Upon  combining  these  several  results,  we  obtain 

Vi*)  ({£-£,  M)  /?(£*,£)  .  (42) 

(We  note  that  U(£ts~  Rfo)  '  </(*)  .)  We  see  that  each  beam  ion,  prior  to 
an  elastic  collision,  traverses  the  target  material  with  its  direction  of 
motion  unaltered  but  with  its  energy  diminished  by  electronic  losses 
according  to  the  prescription  Es  £ ’*(*)  .  Elastic  collisions  lead  to  the 
attentuation  of  the  *K  component  of  the  flux  by  the  factor  £(£»,&)  . 

Finally,  we  note  that  we  can  set  £**£?  with  no  loss  of  generality. 


V.  THE  SOLUTION  FOR  ?,(£,*,*) 

If  we  substitute  the  result  (42)  into  Eq.  (25),  we  obtain 

W.  =  ><  fO- <**)  R(E*,e,to)  Ttew  s)  m  (4 

We  can  now  solve  Eq.  (23)  for  by  using  the  method  of  characteristics . 
The  procedure  is  virtually  identical  to  that  used  for  the  case,  so  that 
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we  obtain 


V (£.',*)  =  J“(e)  j**  A  /?(*,£)  Us.  (t, *(<),,# J 

(44) 

where 

X(tJ  =  x  -/<. 

(45) 

XM  -  x  Lt*-e*)/K  . 

(45') 

The  function  t,  =  €»(X(<J)  ,  which  replaces  £»(*)  in  W0  in 
of  Eq.  (44),  is  determined  by 

the  integrand 

X(fJ  -  f.  (€,) 

(46) 

S,  =  [  (f* -£*)]* 

(46') 

where  C  and  &  are  given. 

The  integrand  in  Eq.  (44)  contains  the  factor 
the  expression  (11)  for  ,  and  defining 

fJj  .  By  using 

h  -  Ct/c,)*  t 

(47) 

we  find  that  the  values  of  h  that  satisfy  the  delta  function 

are 

^  -  fit  -  O  t  (S +#'-/)*  )/(*+')  . 

(48) 

The  appropriate  sign  to  be  used  in  Eq.  (48)  for  h  is  indicated  in  Table  1; 
note  that  /<  is  double  valued  as  a  function  of  h  for  the  case  .  The 

quantity  h  is  either  positive  or  zero,  and  is  only  zero  for  the  case  of 

90°  scattering  with  M*t.  The  value  T,  of  €,  corresponding  to  ft  is 

obtained  from  Eqs.  (47,48);  that  is. 

x  -  (  f(ej  - 

(49) 

r.  =  £  (ex V*  -x<  £*)/(/-**)]*. 

(49') 
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Table  1 


Dependence  of  h  ony*  and  AJ 


A1 

— 

e 

A 

h 

>1 

BB 

~!  *  /*  4  / 

i  1 

-/ 

%S  6  >  O 

0  *  /*  4  / 

O  £.  h  *yU  4  1 

<  / 

0  £  es  Om 

f  2  A  £  O-M*)* 

£/,*A.  £  fi-t 

0m  *0  2  0 

(/-#*)*  £A  H 

naam 

(  9m  -  (/-/f3)*  )  (ht  ~  (/<* ) 


We  can  now  write 

U/» t))  ~  P(E,*s/*)  s(t-t)  (50) 

--  0-l(i-s‘t*r'<*)SGO)(i‘~  %r)\  ' ,  (5„ 

where 

C  ~  J*  £,  ,  (52) 

By  combining  these  results  we  obtain 

tCe,,, A)  =/*  S-'(e)  u(ft-rj  <J(e, -7.)  (AT,)  ui^'U-e))  (53) 

'  R(e,,  ij  X(l,  z)  R(i,e)  P(e,r,s<) 

where  *t  and  i  are  given  by  Eqs.  (49,52)  and  where 


13 


(54) 


r 


X  -  X  -/*  ( f  (?)  ~  f(*)) 

X  * 

The  physical  interpretation  of  the  expression  for  is 

straightforward.  Suppose  we  ask  what  the  once  scattered  flux  is  at  depth  x 
with  energy  f  and  direction  cosine  //.  Given  a  particular  form  for  S(£) , 
there  is,  at  most,  one  way  that  a  beam  particle  can  scatter  just  once  and 
have  the  coordinates  (?,*//*)  .  Specifically,  it  travels  without  scattering 
from  to  >)  and  there  scatters  to  after  which  it 

travels  without  further  elastic  collisions  to  The  step 

function  l/C*"' (x-U)  ,  which  arises  from  Eq.  (54),  reflects  the  fact 
that  i  must  be  greater  than  or  equal  to  £.  The  step  function  U(i.)  ,  which 
arises  from  Eq.  (49),  merely  indicates  that  the  beam  particle  must  not  have 
come  to  rest  Defore  it  scattered.  The  explicit  presence  of  these  step 
functions  is  useful  when  evaluating  integrals  over  <*). 


VI .  THE  SOLUTIONS  FOR  %  . . .  ,  1„ 

The  formal  solutions  for  n-l  are  the  same  as  the  first 

steps  in  the  solution  for  #  .  The  method  of  characteristics  yields 

8  (‘.’,+1  =  S  -W  fj'/t  KU,  e)  m^(t,  XM,  />)  ( 55 ) 

where  X(t)  is  again  given  by  Eq.  (45).  For  these  cases  the 
functions  contain  no  delta  functions,  as  W~t  and  ^  did,  so  that  each 
must  be  evaluated  numerically. 

Similarly,  the  expression  for  £&W  is 

i.(e,*,s)  -  s-'(e)  f** (56) 

which  is  an  integral  equation  for  X*. 

Generally,  one  would  not  continue  the  iteration  process  to  large  n  in 
order  to  solve  for  iP  deep  in  the  target;  too  much  computer  time  would  be 


needed  to  calculate  the  %  quantities.  The  main  purpose  of  the  expansion  is 
to  evaluate  the  first  few  terms,  so  that  the  remaining  portion  of  the 
distribution  function,  which  is  found  be  solving  an  integral  equation,  no 
longer  possesses  the  singular  behavior  of  the  first  few  terms. 


VI 1.  SAMPLE  CALCULATIONS 

We  have  evaluated  the  'K  and  contributions  to  $  for  a  number  of 
values  of  the  beam  energy  and  of  the  mass  ratio.  We  will  present  the 
results  for  one  case  in  some  detail  and  then  take  note  of  the  differences 
that  arise  when  M  and  £$  are  varied.  The  purpose  of  these  calculations  is 
to  illustrate  the  sharply  peaked  behavior  of  ¥  near  the  target  surface. 

For  the  elastic  cross  section  we  use  the  form  due  to  Lindhard,  Nielsen 
15 

and  Scharff  (LNS),  which  is  discussed  in  Appendix  B.  For  the  inelastic 
losses  we  use  the  velocity  stopping  law  of  Eq.  (28).  The  constant  K  in  this 
equation  depends  on  the  charges  and  masses  of  the  interacting  ions.  The  LSS 
expression  for  K-,  as  well  as  their  criterion  for  the  applicability  of  this 
law,  are  also  given  in  Apprendix  B.  The  relations  between  if,  7,  and  7 
and  £;  /,  and  /u  can  be  found  in  Eqs.  (48,  49',  52,  54')  with  £,  Or)  given  by 

0,(x)  *  .  (57) 

Unless  otherwise  noted,  all  energies  are  in  units  of  MeV  and  all  depths  in 
microns . 

We  consider  the  case  of  100  keV  helium  ions  incident  on  a  pure  nickel 
target.  We  take  the  energy  transfer  cutoff  77,  which  is  discussed  in 
Appendix  8,  to  be  25  ev. 

Figure  la  shows  the  region  of  the  plane  for  which  is 

non-zero;  the  restrictions  that  form  the  boundaries  in  the  figure  are  also 
shown.  These  boundaries  arise  from  the  step  functions  that  appear  in 
Eq.  (53).  The  plot  shown  is  for  a  depth  of  0.04  microns  (£>  =  87.998  keV), 
but  the  results  at  other  depths  are  similar.  Figure  lb  is  a  magnified  view 
of  the  small  angle  portion  of  Fig.  la.  We  see  that  the  77  cutoff  allows  a 
minimum  scattering  angle  of  9mtm-  3.47°.  More  generally,  this  minimum 
angle  is  given  by 
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(a)  Region  of  the  plane  for  which  Wfw)  is  non-zero  for 
iOO  keV  He  incident  on  Ni  at  X  =  .04.  (b)  Expanded  view  of  the 

small  angle  portion  of  (a). 


Om.„  =  ~'L  i  0¥+f)O-  T, /*J  -  *(#’0  (/-  T,/f,)  J 


(58) 


We  next  present  a  plot  of  the  logarithm  of  the  absolute  value  of  the 
quantity 

//  (*,*)  “  ff°  Wfe**At)  (59) 

which  measures  the  particle  flux  of  all  energies  at  the  depth  x  and  in  the 
direction  QsfS'/A.  .  Figure  2  shows  that  almost  all  of  the  Vf  flux  is 
concentrated  at  very  small  angles  and,  consequently,  the  flux  is 
concentrated  in  a  very  narrow  energy  range.  In  Fig.  3,  we  show  an  isometric 
plot  of  the  logarithm  of  for  CTi  9*f04  at  jr  =  0.04  microns.  We 

again  see  the  strong  peaking  for  small  scattering  angle.  In  the  region 
where  £  can  approach  zero  (  blows  up  as  the  minus  one  half 

power  of  £,  which  is  depicted  by  the  arrow  rising  from  this  region  in 
Fig.  3.  The  constant  of  proportional ity  is  so  small,  however,  that  we  can 
infer  that  only  a  negligible  faction  of  the  total  flux  comes  to  rest  at  this 
depth. 

We  have  made  similar  calculations  of  for  a  variety  of  values  of  the 
mass  ratio  and  of  the  beam  energy  £».  The  results  are  as  expected.  As  the 
beam  energy  increases,  the  peaking  also  increases,  due  in  part  to  the 
decrease  in  the  minimum  scattering  angle.  As  the  incident  particle  becomes 
more  massive,  in  comparison  to  the  target  particles,  the  flux  again  is  more 
strongly  peaked  in  the  forward  direction.  This  is  due  both  to  the  reduction 
in  the  allowed  minimum  scattering  angle  and  to  the  fact  that  the  kinematics 
favor  smaller  angle  scattering  as  ft  decreases. 

In  general,  the  single  elastic  scattering  contributions  to  $(2,^/*)  is 
sharply  peaked  in  the  forward  direction,  with  most  of  the  flux  concentrated 
at  angles  near  the  minimum  scattering  angle.  This  angle  is  determined  by 
the  cutoff,  which  precludes  elastic  collisions  in  which  the  primary  knock-on 
atoms  would  have  less  energy  than  is  necessary  to  allow  them  to  escape  from 
their  lattice  sites. 
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180  150  120  90  60  30  0 


9  =  COS-1  jJL 

Fig.  2.  Plot  of  log  \l>  fr/W/  vs.  angle  for  100  KeV  He  incident  on  Ni 
at  x=0.04. 
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9  =  cos 1  /x 

Fig.  3  Plot  of  log  vs.  angle  for  100  keV  He  incident  on  Ni 

at  x«0.04. 
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Vlll.  CONCLUSIONS 


The  present  work  can  serve  as  the  starting  point  for  a  number  of 

investigations.  One  could  solve  for  by  iteration,  that  is,  by 

calculating  successive  ^  terms.  Near  the  surface,  only  a  few  iterations 

would  be  needed,  so  that,  for  instance,  one  could  obtain  a  reliable  estimate 

of  the  distribution  of  beam  ions  that  are  scattered  out  of  the  target.  Some 

account  of  true  surface  effects  would  be  needed.  It  would  also  be 
5 

straightforward  to  include  the  distribution  function  for  the  target 
atoms.  This  distribution  would  have  a  term,  analogous  to  our  term,  that 
described  exactly  the  distribution  of  target  atoms  that  had  been  struck  by 
beam  ions,  but  had  suffered  no  further  elastic  collisions  with  other  target 
atoms.  Successive  contributions  to  the  target  atom  flux  could  be  found  by 
iteraction.  Near  the  surface,  this  calculation  would  yield  the  distribution 
of  sputtered  atoms.  In  principle,  one  could  solve  for  ¥  at  all  depths  by 
iteration,  but  the  cost  in  computer  time  and  the  danger  of  accumulating 
errors  makes  this  approach  unworthy  of  the  effort. 

The  more  common  approach  would  be  to  expand  Im  in  series  of  orthogonal 
polynomials;  Legendre  polynomials  in  /<  and  Hermite  polynomials  in  X,  for 
example.  This  method  would  replace  the  integral  equation  (56)  with  a 
coupled  set  of  much  simpler  integral  equations  for  the  coefficients  of  the 
expansion.  Another,  related,  approach  may  prove  valuable.  We  noted  in  the 
last  section  that  the  %  contribution  to  the  flux  is  strongly  peaked  in  the 
forward  direction.  Clearly,  the  contribution  will  also  be  peaked  in  the 
forward  direction  for  energies  near  the  beam  energy.  Given  this  situation, 
it  may  be  useful  to  expand  as  a  function  of  in  a  series  of  Jacobi 

polynomials.  These  polynomials,  which  are  a  generalization  of  the  Legendre 
polynomials,  have  the  property  that  an  inherent  peaking  can  be  built  into 
them.1*’  By  incorporating  the  peaking  behavior  directly  into  the  expansion 
functions,  we  would  hope  to  reduce  considerably  the  number  of  terms 
necessary  to  represent  the  distribution  function  accurately. 
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APPENDIX  A.  CONTINUOUS  ENERGY  LOSS  MOOEL  FOR  INELASTIC  COLLISIONS 

From  Eq.  (7)  and  the  discussion  following  it  we  see  that  the  inelastic 
(.electronic)  cross  section  furnishes  two  contributions  to  the  Boltzmann 
equation;  one  arises  from  the  total  cross  section  term  on  the  left-hand  side 
of  the  Boltzmann  equation  ano  the  other  from  the  scattering  term  on  the 


right-hand  side.  We  write  these  as 

Cl  *  (Al) 

Cfi  -  *7*  (£'■**)  (./•  W ''  ¥(£,*,/*')  (A2) 

where 

(f'i  s  Jd*!/  97*  (A3) 

As  indicated  in  Eq.  (10),  we  wish  the  inelastic  cross  section  to  satisfy 

s(6')  -  (.£'•£)  (A4 ) 

wnere  5(0)  is  the  rate  of  energy  loss  of  beam  ions  to  purely  inelastic 

...  .  10 
collisions . 

Following  Symon^  we  use  a  continuous  energy  loss  model.  A  general 


form  for  vj,  (  a'**)  which  can  be  made  to  satisfy  Eqs.  (A3,A4)  is 

lilt'**)  =  TZ<m(0')  %•''*')  +f'(r)  jyUw)~]  (as) 

where  /  and  are  yet  to  be  determined.  Upon  substituting  the  form  (A5) 
into  Eq.  (A3)  we  find 

/  -  l  6(0)  SC*'*')  *  6(0') £  £(*-*')  J  #  (A6 ) 

We  can  evaluate  the  integral  in  Eq.  (A6)  by  using  the  relations 

X#  **  i(x)  H*-y)  -  ftyJ  (A7 ) 
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/-  A  j(*)  £  S(/-y)  =  *<*'/)  j'(x)  =  -j'ffJ  (A8) 

so  that  we  obtain 

/  •  r*  4(e'J  m 

Thus  ’&(&!  is  determined  as  a  functional  of  /,(ej  ,  which  is  still 
undetermined . 

After  substituting  the  form  (A5)  into  Eq.  (A4),  we  obtain 

Ste'j  -  M »,  v3  77.{0(e'J  (aio) 


In  similar  fashion  we  can  evaluate  obtaining 

Cfi  -  At[  v  £  (  V 

*  v*  r, rtm(g)  S>(£)  ?(£,*,*)  ] 


(All) 


which  becomes,  upon  inverting  the  result  (A9), 

-  3/s  *  **.,'•(£)  /,(£)  ¥ (?.***) 

,  /  v  ,  ( A 1 2 ) 

+  Vv*  5*  (  %tn(S)  /,(£)  $(&,*))  +  A/  J5t /,/&  ?(&/*). 

By  combining  Eq.  (AIO)  with  Eq.  (A12)  we  find  that 

cn  =■  v  r,tim(e)  •Jte./s')  +  h  (  (Ai3) 

We  see  that  when  and  4*  are  substituted  into  the  Boltzmann  equation,  the 
first  term  in  cancels  C£,  so  that  the  final  form  for  Boltzmann  equation 
is  that  given  in  Eq.  (11). 

We  note  that  *£,,*,  which  in  the  case  of  elastic  scattering  would 
correspond  to  the  total  cross  section,  is  an  arbitrary  quantity,  depending 
on  the  choice  of  the  function  /,(&).  We  conclude  that  the  concept  of  a 


23 


1 


total  cross  section  has  no  meaning  for  a  differential  cross  section  of  the 
form  (A5). 


APPENDIX  B.  LINDHARD  CROSS  SECTION 

15 

The  lindhard,  Nielsen,  and  Scharff  (LNS)  differential  cross  section 
for  elastic  atomic  scattering  is  given  by 

AbnW-v)  S« *J I -W'l o»r‘ j* st  ( 

where 

t  *  e' ( e)  /  ( r  A* )  . 

The  Lindhard  unit  of  energy  is  given  by 

*2,2*  (tYl) 

*  0. 1fS3  4.  3~* 

2  =  (g,*  -  M**)* 

/  -  t",  *fL  (m,  +  ) 

where  is  the  Bohr  radius  for  the  hydrogen  atom.  The  function  /«*)  Is 

obtained  by  LNS  from  a  Thomas-Fermi  model  of  atomic  scattering  and  is 
supplied  by  them  in  tabular  form. 

When  the  collision  kernal  f(e:e)  is  expressed  in  terms  of  the  LNS  cross 
section  we  obtain 

F(e'fe)  r  i  rr<  i/'J  #t*)  \*/*\  ( b7 ) 

so  that 

~  (m)’*  S(.t  F(f,e)  Sr  .  (B8) 


( B2 ) 

(B3) 

(B4) 

(B5) 

(B6) 
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The  total  cross  section  is  given  by 


Tf(e')  -  v  +*  Jf/  S?  ?’J  /(f)  (B9) 

where 

t,  *  (e'v/0‘f.'))t  ,  fj  -  e'/fi.  .  (bio) 

Since  the  LNS  cross  section  blows  up  as  the  energy  transfer 
goes  to  zero,  it  is  necessary  to  introduce  a  minimum  energy  transfer  7^. 

The  imposition  of  the  cutoff  T,  is  less  arbitrary  then  it  might  seem  since 
7/  can  be  related  to  the  displacement  energy  which  is  the  average  energy 
needed  to  displace  a  target  atom  from  its  lattice  site. 

The  function  /w/ is  presented  in  tabular  form  by  LNS,  but  it  is  more 
convenient  to  use  the  analytical  fit  to  ft,*)  provided  by  Winterbon, 

Sigmund,  and  Sanders  (WSS).  Specifically,  they  give 

/(■**>  =  3  -f  *[/r  ax  **)*]'*  (Bn) 

with  a**?*T.  An  explicit  expression  for  the  total  cross  section 
corresponding  to  the  WSS  choice  for  /«*)  has  been  developed.^7 
LSS1^  provide  an  expression  for  the  constant  that  appears  in  the 
velocity  stopping  law.  If  we  write 

F(Ct  5  t<  F *  ( B 1 2 ) 

with  $(£)  in  units  of  MeV  per  micron  and  F  in  MeV,  we  have 

J<  -  /«?/*  ■  /°~33  2,*  (2,27/2)  A/ ((*’*)  ( («*«))  ^  .  (B13) 


LSS  give 

B  <  6.0JS"  2,^  (*•*«)  (B14) 

as  a  rough  criterion  for  the  applicability  of  the  stopping  law  ( B1 2 ) . 
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